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The hadronic decays of non-strange negative parity baryons are analyzed in the framework of 
the 1/N C expansion. A complete basis of spin-flavor operators for the pseudoscalar meson S and 
D partial wave amplitudes is established to order 1/N C and the unknown effective coefficients are 
determined by fitting to the empirically known widths. A set of relations between widths that result 
at the leading order, i.e. order N®, is given and tested with the available data. 
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IT) 1 I- INTRODUCTION 

o : 

The 1/N C expansion has proven to be a very useful tool for analyzing the baryon sector. This success is mostly a 
consequence of the emergent contracted spin-flavor symmetry in the large N c limit [1, 2]. For ground state baryons 
(identified for jV c = 3 with the spin 1/2 octet and spin 3/2 decuplet in the case of three flavors), that symmetry 
gives rise to several important relations that hold at different orders in the 1/N C expansion [3-5]. The domain of 
excited baryons (baryon resonances) has also been explored in the framework of the 1/N C expansion[6]~ [f 4] with very 
promising results. The analyses carried out so far have been constrained to states that belong to a definite spin-flavor 
and orbital multiplet, i.e. the possibility of mixing of different such multiplets (so called configuration mixing) has 
been disregarded. It is very likely that such effects are small for dynamical reasons. Indeed, it has been shown that 
the only configuration mixings that are not suppressed by 1/N C factors involve couplings to the orbital degrees of 
freedom[15]. The 1/N C analyses of excited baryon masses have shown that orbital angular momentum couplings turn 
out to be very small [8, 9], which is in agreement with older results in the quark model[16]. This strongly suggests that 
t-H ' a similar suppression, which is not a consequence of the 1/N C expansion but rather of QCD dynamics, also takes place 
in configuration mixings. Thus, disregarding configuration mixing is likely to be a good approximation for the purpose 
of phenomenology. Within such a framework, a few analyses of excited baryon strong decays have been carried out, 
namely the decays of the negative parity SU(6) 70-plet[10] and of the Roper 56-plet[lf]. In the case of interest in the 
present work, namely the 70-plet, the analysis in Ref. [10] used an incomplete basis of operators at sub-leading order 
' in 1/N C . The main motivation of the present work (see Ref. [17] for more details) is to provide a general framework 
for the study of the strong decays and to perform a complete analysis to 0(1/ N c ) in the particular case of the decays 
O ' of the non-strange members of the 70-plet (i.e., the mixed symmetry 20-plet of SU(4)) into ground state baryons plus 
D , a pion or an eta meson. 

' This contribution is organized as follows: Section 2 contains the framework for calculating the decays, Section 
3 provides the basis of effective operators, Section 4 presents the results, and finally the conclusions are given in 
Section 5. 
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II. FRAMEWORK FOR DECAYS 



As discussed above in the present application of the 1/N C expansion to excited baryons the assumption is made 
that there is an approximate spin-flavor symmetry. Thus, the excited baryons are therefore classified in multiplets 
of the 0(3) x SU(2Nf) group. 0(3) corresponds to spatial rotations and SU(2Nj) is the spin-flavor group where 
Nf is the number of flavors being considered, equal to two in the present work. The ground state baryons, namely 
the N and A states, belong to the (1, 20s) representation, where the 20s is the totally symmetric representation of 
SU(4). The negative parity baryons considered here belong instead to the (3,20ms) representation, where 20ms is 
the mixed symmetric representation of SU(4). For general N c the spin-flavor representations involve, in the Young 
tableaux language, iV c boxes and are identified with the totally symmetric and the mixed-symmetric representations 
of type (N c — 1, 1) for ground and excited negative parity states respectively. Since in the mixed symmetric spin- flavor 
representation one box of the Young tableaux is distinguished, such a box is associated with the "excited quark" in 
the baryon. In a similar fashion, and without any loss of generality it is possible to distinguish one box in the ground 
state multiplet as well. This is a very convenient procedure that has been used repeatedly in previous works. The 
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spin and isospin quantum numbers of the distinguished box will be denoted with lower cases, and the corresponding 
quantum numbers of the rest of N c — 1 boxes (which are in a totally symmetric representation of SU (4) and form the 
so called "core" of the large iV c baryon) after they are coupled to eigenstates of spin and isospin will be denoted by 
S c and I c respectively. Notice that S c = I c for totally symmetric representations of SU(A). For a given core state, 
the coupling of the excited quark gives eigenstates of spin S and isospin I: | S, S3; 1 , I 3 ; S c ) . These states are not 
in an irreducible representation of SU(4), as they arc not in an irreducible representation of the permutation group. 
The totally symmetric states are given by: 

\S,S 3 ;I,I 3 ) S = Yl C S( S >V) |S,S 3 ;/,/ 3 ;S c = S + 77), (1) 

r,=±± 

where 



r ,o+S / (25+l±l)(iV c + lT(2S + l)) 

C S {S,±j = J 2N c (2S+l) ' (2) 



while the mixed symmetric (MS) states (N c — 1, 1) are given by: 



S,S 3 ,II 3 ) MS = C MS (I,S, V )\ S,S 3 ;I,I 3 ;S c = S + ri), (3) 

r,=±i 



where 



f 1 if I = S±1 

tWI.S.iJ) = OiiI = S T l (4) 

2 , / (2S+l T l)(jV c + l±(25+l)) - f r _ a 

y 2N C (2S+1) 

Finally, upon coupling the orbital degrees of freedom, the excited baryons in the (3, MS) representation, 
I (1, S) J, J 3 ; I, I3; S)-^c, are obtained. 

Note that there are two sets of N* states each consisting of two states with the same spin and isospin. The physical 
states are admixtures of such states, and are given by: 



N} \ _ ( cos 9 2 j sme 2J \ ( 2 NJ 



N*j ) \ -sin6> 2 j cos6» 2 .7 J \ 4 iV, 



4 7\r* 



(5) 



1 3 *( f ) 

where J = ^ and |, Nj ' are mass eigenstates, and the two mixing angles can be constrained to be in the interval 
[0, 7r). Here the notation 2S+1 7V} has been used. 

Using the standard definition for the decay width and averaging over the initial- and summing over the correspond- 
ing final-baryon spins and isospins, the decay width for each orbital angular momentum lp and isospin Ip of the 
pseudoscalar meson is given by 



r [/p,Jp] = j ps 2^°q n q yip,i P ,D,i,j ,i ,j )\ ^ 



E 9 Cj tp ' lp] B q (£p,I P ,S,I,J* , I* , S* ) 1 2 
V'(2J* + 1)(2J» + 1) 

where the phase space factor f ps is 



f -JlP ( 7 ) 

Here, we have used that the baryonic operator admits an expansion in 1/N C and has the general form: 



B i£pJp] _ 



where 

K:«f Pl ) = e^^j.j* 1 & . ( 9 ) 
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and the factor (^f ) P is included to take into account the chief meson momentum dependence of the partial wave. 
The scale A is chosen in what follows to be 200 MeV. Here, is an operator that produces a transition from the 

triplet to the singlet 0(3) state, and (^Gy' Ip ^j is a spin-flavor operator that produces the transition from the mixed- 
symmetric to the symmetric SU (4) representation. The label j denotes the spin of the spin-flavor operator, and as 
it is clear, its isospin coincides with the isospin of the emitted meson. The dynamics of the decay is encoded in the 
effective dimensionless coefficients cf F ' Ip ^ {kp). The reduced matrix elements B q (£p, Ip, S, I, J* ,1* , S*) appearing in 
Eq.(6) can be easily calculated in terms of the reduced matrix elements of the spin-flavor operators. Note that in the 
present case, where ip can be or 2 only, the spin- flavor operators can carry spin j that can be 1, 2 or 3. 

The terms in the right hand side of Eq.(8) are ordered in powers of 1/N C . As it has been explained in earlier 
publications [3], the order in l/N c is determined by the spin-flavor operator. For an n-body operator, this order is 
given by 

v = n — 1 — k, (10) 

where k is equal to zero for incoherent operators and can be equal to one or even larger for coherent operators. More 
details can be found in the following section where a basis of operators G is explicitly built. 

With the definition of effective operators used in this work, all coefficients Cq' F,Ip]> (kp) in Eq.(8) are of zeroth order 
in N c . The leading order of the decay amplitude is in fact N® [15]. At this point it is important to comment on the 
momentum dependence of the coefficients. The spin-flavor breakings in the masses, of both excited and ground state 
baryons, give rise to different values of the momenta kp. In this work, we have adopted a scheme where the only 
momentum dependence assigned to the coefficients is the explicitly shown factor (kp/A) lp that takes into account 
the chief momentum dependence of the corresponding partial wave, and the rest of the dependence is then encoded 
in the coefficients of the sub- leading operators. 



III. BASIS OF SPIN-FLAVOR OPERATORS 



The spin-flavor operators to be considered in the present paper must connect a mixed-symmetric with a symmetric 
representation. Generators of the spin-flavor group acting on the states obviously do not produce such connection. 
However, generators restricted to act on the excited quark or on the core of iV c — 1 quarks can do this. The spin flavor 
operators can, therefore, be represented by products of generators of the spin-flavor group restricted to act either on 
the excited or on the core states. In the following the generators acting on the core are denoted by S c , G c , T c , and 
the ones acting on the excited quark by s, g, t. The generators G c are known to be coherent operators, while all the 
rest are incoherent. In order to build a basis of operators for the present problem one has to consider products of 
such generators with the appropriate couplings of spins and isospins. The n-bodyness (nB) of an operator is given by 
the number of such factors, and the level of coherence of the operator is determined by how many factors G c appear 
in the product. It should be noticed that in the physical case where N c = 3, only operators of at most 3B have to 
be considered. Still, to order 1/N C there is a rather long list of operators of given spin j and isospin Ip. This list 
can be drastically shortened by applying several reduction rules. The first rule is that the product of two or more 
generators acting on the excited quark can always be reduced to the identity operator or to a linear combination 
of such generators. The second set of rules can be easily derived for products of operators whose matrix elements 
are taken between a mixed-symmetric and a symmetric representation. These reduction rules are as follows (here A 
represents generators acting on the excited quark and A c represent generators acting on the core): 



A - -A c 

(A c )i(A c ) 2 = -Ai(A c ) 2 - A 2 (A c )i + IB operators. (11) 
Therefore, only the following types of operators should be considered 

IB 2B 3B 

A ; ^ A i(A c ) 2 ; ^ Ai (A c ) 2 (A c ) 3 (12) 

It is convenient to make explicit the transformation properties of each basic operator under spin j and isospin t. In 
what follows we use the notation 0b'.*] to indicate that the operator O has spin j and isospin t. It is easy to see that 

X M = s [i,o ]; t [o,i] ) ff [i,i] 

(A c )^ = (S C )W, (T C )[M 5 (Gc) [i,i] (13) 
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For decays in the 77— channels the spin-flavor operators transform as \j, 0] while for decays in the pion channels they 
should transform as [j, 1], where in both cases j = 1,2,3. Knowing the transformation properties of each basic operator 
given in Eq.(13) it is easy to construct products of the forms given in Eqs.(12) with the desired spin and isospin. 
Thus, the IB operators that contribute to a given [j,t] are just those A^'*! given in Eq.(13) which have the proper 
spin and isospin quantum numbers. Similarly the possible 2B operators are given by the products A^ 1 '* 1 ' (A c )!f '* 2] 
coupled to the required [j, t] by means of the conventional Clebsch-Gordan coefficients. 

In order to construct all possible 3B operators, it is convenient to note that there are additional reduction rules 
for the products (A C ) 2 (A C ) 3 . The starting point is to consider all possible products of two core operators. Since 
one is interested in keeping contributions of at most order 1/N C , and because these products will appear only in 3B 
operators, at least one of the operators must be a G c . Using the reduction relations[2], it can be shown[17] that the 
relevant list of independent products of two core operators turns out to be 

({T C ,G C })^ , {[S c ,G c ]fX , ({S c ,G c }f^ , (G c G c )™ (14) 

where ([T c , Cc])' 1 ' 1 ' denotes (TcGe)! 1 ' 1 ' — (GcTc)^ 1 ' 1 ^, etc.. By coupling any of these operators with one of the excited 

core operators A^' 1 '* 1 ' (see Eq.(13)) to the required [j,t] all the possible 3B operators are obtained. 

Using this scheme to couple products of generators it is straightforward to construct lists of operators with spin 
1, 2 and 3 and isospin and 1. Further reductions result from the fact that not all the resulting operators are 
linearly independent up to order 1/N C . The determination of the final set of independent operators for each particular 
decay channel is more laborious since it requires the explicitly calculation of all the relevant matrix elements. The 
resulting basis of independent operators (o^^j^ is shown in Table 1, where for simplicity the corresponding spin 

and isospin projections have been omitted. Finally, the normalized basis operators (B^pil ]j arc defined. They 

differ from those listed in Table 1 by a normalization constant which is determined by requiring that, for N c = 3, 
their largest reduced matrix element should be equal to one for order N® operators and equal to 1/3 for order 1/N C 
operators. 



TABLE I: Operator basis 



n-bodyncss Name Operator Order in 1/N C 

IB op 1 (iW^ 

Pion W ] £ (€(.rc)^)!^ "T- ~ 

Swave 2B £ (( it S c )™) ^ 

olo.i] 1 t c ,„ o ^1,11 V - 1 ' 



-,[2,1] 



s (s(^) [mi )^: 



1 



or ] ifc(f(**) [Mi )K 

Dwave 2B O™ £ (* ( 9 S C )M) ™ 1 



% , a\ 





3B 








({ScGc})^) 12 ' 11 )" 1 


1 












({flcGc})^) 18 - 11 )^ 

/ \i,a] 


1 


Eta 


IB 


O l0,0] 











S wave 


2B 


O [0,0] 




1 

JVc 


(€(.*)^)K:3 


1 


Eta 


IB 


O l2,0] 











D wave 


2B 


n [2,0] 
^2 




1 


(«<-*> [Mi c 


1 






O [2,0] 
^3 




1 

N c 




1 
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IV. RESULTS 



The different empirical S- and D-wave decay widths used in the analysis are the ones provided by the PDG[18]. 
The values for the widths and branching ratios are taken as the ones indicated there as "our estimate" , while the 
errors are determined from the corresponding ranges. The corresponding partial widths obtained from those values 
are explicitly displayed in Table 3. The entries indicated as unknown reflect channels for which no width is provided 
by the PDG or where the authors consider that the input is unreliable, such as in the 7rA decay modes of the iV(1700) 
and 7V(1675) and the D-wave r\— decay modes. At this point it is important to stress the marginal precision of the 
data for the purposes of this work. This work performs an analysis at order 1/N C , which means that the theoretical 
error is order 1/N%- This implies that amplitudes are affected by a theoretical uncertainty at the level of 10%. Thus, 
in order to pin down the coefficients of the subleading operators, the widths provided by the data should not be 
affected by errors larger than about 20%. As shown in Table 3, the experimental errors are in most entries 30% or 
larger. In consequence, the determination of the subleading effective coefficients is affected by large errors as the 
results below show. 

Before presenting the results of the fits, it is convenient to derive some parameter independent relations that can 
be obtained to leading order. These relations serve as a test of the leading order approximation. Since at this order 
there are only four coefficients and two angles to be fitted, and there are a total of twenty partial widths (excluding all 
D-wave rj channels but including kincmatically forbidden ^-channel decays) , there are fourteen independent parameter 
free relations that can be derived. These relations are more conveniently written in terms of reduced widths, i.e., 
widths where the phase space factor f sp (see Eq.(7)) has been removed and denoted here by V. Considering the S-wave 
decays in the n mode, there are six decays and three parameters in the fit. Thus, three parameter free relations must 
follow. These relations and the corresponding comparison with experimental values read: 



^(1535) — ttN . r A r( 152 o)^ 7r A 



A(1620)-^7rJV : J- A(1700)^ttA 

(15) 



+rjv(i65o)->7riv +rjv(i 700) — »7rA 

Th. 1:1: 0.17 : 0.42 

Exp. 1 : unknown : 0.19 ±0.07 : 0.62 ± 0.33 



Within the experimental errors the relations are satisfied. In a similar fashion, relations involving D- wave decays in 
the 7r mode can be obtained. There are in this case four parameters to fit and eleven partial widths. Thus, seven 
parameter free relations can be obtained. In general these relations are quadratic and/or involve some of the unknown 
decay widths. However, three testable linear relations can be obtained. The simplest of them reads, 



2 rAflfi2(VI-nrA + rAfi7nni->irA — 8T 
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A(1620)^ttA T -1- A(1700)^ttA — O 1 A(1700)^7rAf T — 1 jV(1675)^7rJV 

Exp. 5.24 ±1.95 = 2.19 ±0.62 

As we see this relation is not well satisfied by the empirical data. Something similar happens with the other two. 
Thus, one can anticipate a poor leading order fit to the D-wave pion decays. 

In Table 2 the results of several fits are displayed. In these fits the decay amplitudes are expanded keeping only the 
terms that correspond to the order in 1/N C of the fit. Similarly, when performing the fits the errors have been taken 
to be equal or larger than the expected accuracy of the fit (30% to the LO fits and about 10% for the NLO ones). 

The first LO fit only considers the S-wave 7r-modes. One notices here that both 6\ and 63 have a two fold ambiguity 
at this order. The second leading order fit includes the D-waves and 77-modes. The angles remain within errors equal 
to the ones from the first fit. Table 3 shows that the iV(1535) — > r/N width results to be a factor four smaller than 
the empirical one (this having, however, a rather generous error). In the D-waves, several widths involving decays 
with a A in the final state are also too small. The S-wave 7r-modcs are well fitted and there is no real need for NLO 
improvement. In the D-wave decays in the 7r channel there are two leading order operators that contribute. The 
fit shows that the IB operator has a coefficient whose magnitude is a factor two to three larger than that of the 
coefficient of the 2B operator. The IB D-wave operator as well as the IB S-wave operator oj ' 1 ' stem from the IB 
coupling of the pion via the axial current. Such a coupling naturally occurs as the dominant coupling in the chiral 
quark model[19]. 

The NLO fit involves a rather large number of effective constants. Since there are only sixteen data available, some 
operators must be discarded for the purpose of the fit. It is reasonable to choose to neglect the 3B operators and the 
subleading S-wave operator for 77 emission. The NLO fit has been carried out by demanding that the LO coefficients 
are not vastly different from their values obtained in the LO fit. This demand is reasonable if the assumption is made 
that the 1/N C expansion makes sense. The fact that the NLO coefficients do not have unnaturally large values with 
respect to the scale set by the LO fit indicates the consistency of the assumption. This clearly is no proof, however, 
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TABLE II: Fit parameters. Fit #1: Pion S-waves LO. In this case there is a four-fold ambiguity for the angles {9\, 83} given by 
the two values shown for each angle. Fit #2: Pion S and D-waves, eta S-waves, LO. In this case there is a two-fold ambiguity 
for the angle 81. For the angle 83 there is an almost two- fold ambiguity given by the two values indicated in parenthesis and 
which only differ in the two slightly different values of Cg 2,1 '. Fit #3: Pion S and D-waves, eta S-waves, NLO, no 3-body 
operators. No degeneracy in 9i but almost two-fold ambiguity in 83 given by the two values indicated in parenthesis. Values of 
coefficients which differ in the corresponding fits are indicated in parenthesis. 



Coefficient 


#1 LO 


#2 LO 


#3 NLO 




31 ±3 


31 ± 3 


23 ±3 




- 


- 


(7.4,32.5) ± (27,41) 








(20.7,26.8) ± (12, 14) 








(-26.3, -66.8) ± (39,65) 






4.6 ±0.5 


3.4 ±0.3 








-4.5 ±2.4 




- 


- 


(-0.01,0.08) ±2 


C P,1] 


- 


- 


5.7 ±4.0 


r o 1 1 

4 1 






3.0 ±2.2 






(-1.86, -2.25) ±0.4 


-1.73 ±0.26 


C [2,l] 

c h,i] 






c [o,oj 




11 ±4 


17 ±4 


c \o,o] 








81 


1.62 ± 0.12 
0.29 ± 0.11 


1.56 ± 0.15 
0.35 ± 0.14 


0.39 ±0.11 




3.01 ± 0.07 
2.44 ±0.06 


(3.00, 2.44) ±0.07 


(2.82,2.38) ±0.11 


^dof 


0.25 


1.5 


0.9 


dof 


2 


10 


3 



TABLE III: Partial widths resulting from the different fits in Table II. Values indicated in parenthesis correspond to the cases 
in which almost degenerate values of 83 lead to different partial widths. 







Decay 




Emp. 
[MeV] 


#1 LO 
[MeV] 


#2 LO 
[MeV] 


#3 NLO 
[MeV] 






7V(1535) -> 


7TN 


68 ± 27 


74 


62 


(58,68) 






JV(1520) -> 


ttA 


10 ±4 


10 


9.7 


9.5 




7T 


Af(1650) -> 


ttN 


121 ± 40 


132 


144 


122 


s- 


wave 


JV(1700) -» 


7rA 


unknown 


175 


175 


(259,156) 






A(1620) -> 


vrAf 


38 ± 13 


35 


35 


38 






A(1700) -> 


7rA 


112 ±53 


81 


81 


(135,112) 






N(1535) -> 


7rA 


1±1 




.01 


0.5 






JV(1520) 


ttAT 


67 ±9 




70 


65 






AT(1520) -v 


7rA 


15 ±3 




2.8 


13 






AT(1650) -> 


7rA 


7±5 




0.12 


8 




7T 


^(1700) -> 


ttAT 


10 ±7 




10 


(11,9) 


D 


-wave 


N(noo) -» 


ttA 


unknown 




4 


(4,9) 






Af(1675) 


vrAT 


72 ± 12 




85 


76 






W(1675) -» 


7rA 


unknown 




45 


79 






A(1620) -> 


7rA 


68 ± 26 




30 


87 






A(1700) -» 


ttAT 


45 ±21 




49 


32 






A(1700) -» 


ttA 


12 ± 10 




15 


18 




V 


AT(1535) -» 


r;JV 


64 ±28 




17 


(57,61) 


S- 


wave 


AT(1650) -> 




11 ±6 




14 


12 
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that the 1/N C expansion is working. As mentioned earlier, the chief limitation here is due to the magnitude of the 
errors in the inputs. This leads to results for the NLO coefficients being affected with rather large uncertainties. 
Indeed, no clear NLO effects can be pinned down, as most NLO coefficients are no more than one standard deviation 
from zero. Because the number of coefficients is approximately equal to the number of inputs, there are important 
correlations between them. For instance, the S-wave NLO coefficients are very correlated with each other and with 
the angle 6*3. For the S- waves the LO fit is already excellent, and therefore nothing significantly new is obtained by 
including the NLO corrections. Correlations are smaller for the D-wave coefficients. As mentioned before, here the 
LO fit has room for improvement, and thus the NLO results are more significant than in the case of the S-waves. 
Still, no clear pattern concerning the NLO corrections is observed. One interesting point, however, is that without 
any significant change in the value of 81 the 77-modes are now well described. The reason for this is that in the LO fit 
the matrix elements of the operator o[ ' ' were taken to zeroth order in 1/N C , while in the NLO fit the \/N c terms 
arc included. These subleading corrections enhance the amplitude for the 2 N^ 2 and suppress the amplitude for the 
4 -^i/2' This along with an increment in the coefficient brings the fit in line with the empirical widths. One important 
point is that at NLO the two fold ambiguity in 8\ that results at LO is eliminated. The smaller mixing angle turns 
out to be selected. The angle #3 remains ambiguous and close to the values obtained in the LO fits. It should be 
noticed that the present values of both mixing angles are somewhat different from the values Q\ — 0.61, 8 3 = 3.04 
obtained in other analyses[10, 16, 20]. 

V. CONCLUSIONS 

In this work we have analyzed the strong decays of the negative parity excited baryons in the 1/N C expansion 
assuming that configuration mixings can be neglected. A basis of effective operators, in which the S- and D-wave 
amplitudes are expanded, was constructed to order 1/N C . All dynamical effects are then encoded in the effective 
coefficients that enter in that expansion. We find that a consistent description of the decays within the l/N c expansion 
is possible. Indeed, up to the relatively poor determination of 1/N C corrections that results from the magnitude of 
the errors in the input widths, these corrections are of natural size. A few clear cut observations can be made. The 
most important one is that the S-wave it- and 77-channels are well described by the leading order operators (one for 
each channel) provided one includes the contributions subleading in \/N c in the matrix elements for the 77-decays. 
The mixing angle Q\ is then determined by these channels up to a twofold ambiguity, which is lifted when all channels 
are analyzed at NLO. The angle 83 is also determined up to a two fold ambiguity at LO. The ambiguity remains 
when the NLO is considered. The subleading operators are shown to be relevant to fine tune the S-wave decays and 
improve the D-wave decays. Because of the rather large error bars in and significant correlations between the resulting 
effective coefficients, no clear conclusions about the physics driving the 1/N C corrections can be made. The mixing 
angles 8 X = 0.39±0.1I and 8 3 = (2.82, 2.38)±0.11 that result at NLO are similar to the ones determined at LO. They 
are , however, somewhat different from the angles 8\ — 0.61 and #3 = 3.04 obtained in other analyses[10, 16, 20]. 
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